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Abstract 



We propose a solar dynamo model distributed in the bulk of the convection 
zone with toroidal magnetic-field flux concentrated in a near-surface layer. We 
show that if the boundary conditions at the top of the dynamo region allow the 
large-scale toroidal magnetic fields to penetrate close to the surface, then the 
modeled butterfly diagram for the toroidal magnetic field in the upper convec- 
tion zone is formed by the sub-surface rotational shear layer. The model is in 
agreement with observed properties of the magnetic solar cycle. 

Subject headings: Dynamo — Magnetohydrodynamics (MHD) — Sun:dynamo 
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Introduction 



It is widely believed that the 11-year sunspot activity is produced and organized by 
large-scale magnetic fields generated somewhere in the deep convection zone. Most of the 
solar dynamo models suggest that the toroidal magnetic field that emerges on the surface 
and forms sunspots is generated near the bottom of the convection zone, in the tachocline 



or just beneath it in a convection overshoot layer, (see, e.g., Choudhuri et al. 1995, Riidiger 



& Brandenburg, 


1995; Dikpati &; Charbonneau, 1999, Bonanno et al. 


2002, 


Tobias & 



Weiss, 2007). The belief in a deep-seated solar dynamo comes from the fact that this region 



is sufficiently stable, and can store magnetic flux despite the magnetic flux-tube buoyancy 



effect (Parker, 1975, Spiegel fc Weiss, 1980, van Ballegooijen, 1982, Spruit fc Roberts, 1983 



van Ballegooijen & Choudhuri, 1988, Choudhuri, 1990). The tachocline represents a strong 



radial shear of the angular velocity. Yet, turbulent diamagnetism (see, e.g., Zeldovich, 1957 
or 



Kitchatinov & Riidiger, 1992) pumps the magnetic fields from the intensively mixed 



interior of convection zone to its boundaries. This effect can substantially amplify the 



toroidal magnetic fields near the convection zone boundaries, (see, e.g., Krivodubskij , 1987 



Guerrero & de Gouveia Dal Pino , 2008 ) 



However, an attention was drawn to a number of theoretical and observational 



problems concerning the deep-seated dynamo models (Brandenburg, 2005, 2006). A 



renewed discussion of the place of the solar dynamo can be found, e.g., in papers by 



Brandenburg (2005) and Tobias & Weiss (2007). In particular, there are some arguments 



that the sub-surface angular velocity shear could play an important role in the dynamo 
distributed in the convection zone. This shear layer becomes an important ingredient of the 



flux-transport dynamo models as well, (see, e.g., Guerrero & de Gouveia Dal Pino, 2008). 



In this Letter we discuss the importance of the surface boundary conditions for the 
dynamo models, which include the subsurface shear layer. The boundary conditions 
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commonly used in the dynamo models correspond to a perfect conductor at the bottom of 
the convection zone and vacuum boundary conditions at the top. Both the vacuum and 
perfect conductor boundary conditions can be regarded as a mathematically convenient 
idealization. The top boundary conditions play a particularly important role because they 
control the escape of the dynamo generated magnetic fields from the Sun. 

The perfect conductor boundary condition is usually identified as "closed", (e.g., 



Choudhuri, 1984), because in this case there is no penetration of the generated magnetic 



flux to the outside. For the axi-symmetric magnetic fields all magnetic field flux is closed 
inside the dynamo region. The vacuum boundary condition is identified as "open". In this 
case the poloidal field lines are open to the outside, and the corresponding poloidal magnetic 
flux "freely" escapes. Also, the strength of the toroidal magnetic field goes smoothly to zero 
at the boundary. This means that the vacuum boundary condition does not allow to the 
toroidal field penetrate to the surface. With such boundary condition it is hardly possible 
to form sunspots from the near-surface large-scale toroidal magnetic fields. 

Bearing in mind the dynamical nature of magnetic fields on the solar surface one 
can model the near-surface behaviour by using a combination of the "open" and "closed" 
types of the boundary conditions. Various consequences of this idea were explored, (see, 
e.g., IChoudhurij [T9841 |Tavakol et alj [T9951 |Covaset al.[ [T998| |Kitchatinov fc Mazur[ [T999 



Kitchatinov et al. 2000, Tavakol et al. 2002, Kapyla et al. , 2010). Here, we apply this 



approach to a solar dynamo model that extends from the bottom of the convection zone 
to the top, including the region of the strong sub-surface rotational shear. We show that 
allowing the toroidal magnetic flux to penetrate to the surface brings the butterfly diagram 
of the toroidal large-scale magnetic field in the upper convection zone and also the phase 
relations between the different components of the dynamo-generated magnetic field in 
agreement with solar-cycle observations. 
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2. Dynamo equations 



The evolution of the axi-symmetric magnetic field (B being the azimuthal component 
of the magnetic field, A is proportional to the azimuthal component of the vector potential) 
is governed by the following equations: 



dA 

~dt 
OB 

~dt 



'dQdA dndA\ 1 dr£ e sin 9 d£ r 
dr djj, djj, dr J r dr r dp ' 



where, 



r sin 9£rk 



= ^ {/«> + 2/W} ^ + ^ V T {ff + 2/«} 
+ iP a C a TtrGf^rfjLsm9f(e)B 
+ G^f^9 d ^-(f^-f^s^9)f) 



[1 - fi 2 ) d 2 A 
r 2 dji 2 



+ 



VTlpr, 



d sin 9B 



dfj, 



2f[ a) a sin 9^-rB + 2rGfF'fJi sin #5 



<9 



+ 



rfc = ^[(/ 2 (d) + 2A (a V 2 
+ 2/ 1 (a) //sin0^-(sin0£) 
- rG (7 1 (a) cos20 + / 3 (a) ) B 



<9r£> 
dr 



(1) 
(2) 



(3) 



(4) 



(5) 



These equations are similar to those used by Pipin & Seehafer (2009) and Seehafer & 



Pipin (2009). We use the same notations for the functions and parameters as in the paper 



of Pipin (2008) (hereafter, P08). Here, G = d r \ogp is the density stratification scale. 
Functions /i^io depend on the Coriolis number fi* = 2r c fi ; functions ip v>a describe 



magnetic quenching and depend on /3 = B/y popu 2 . For reference, these functions are given 
in Appendix. The parameter C a controls the strength of the a-effect. In the presented 
model the a-effect is distributed in the bulk of the convection zone. For a more clear 
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demonstration of the boundary condition impact, we confine the a-affect in a low-latitude 
region where the radial gradient of the angular velocity is positive in the most part of 



the solar convection zone. Similarly to Dikpati et al. (2004) we specify the confinement 
function: 

/ (0) = (l +e 3O(|e- 7r /2|- 7r /6))- 1 . ( 6 ) 

In the radial direction the a-effect depends on the density stratification, G, and function 
of the Coriolis number f$ (O*). We introduce parameter C n to control the turbulent 
diffusion coefficient, rfr = C v r]^\ where rj^ = t c u 2 /3. The internal parameters of the 



solar convection zone are given by Stix (2002). At the top of the solar convection zone 



the stratification is strongly deviates from adiabatic, and also the turbulence parameters 
vary sharply. For this reason we confine the integration domain between O.71_R and 
O.972i? in radius, and it extends from the pole to pole in latitude. The differential rotation 
profile, n = Qofn (x, ji) (shown in Figjlk) is a slightly modified version of the analytical 



approximation proposed by Antia et al. (1998): 

fn(x,fi) = — [fl + 55 (x - 0.7)0 (x,x )<f>(-x, -0.96) 



n 

- 200 (s -0.95)0 {x, 0.96)) 



(7) 



Jp 



+ (2LP 3 (//) + 3P 5 (//)] 



l+exp(^)' 3e 



Jp 



Je [X) 
1 



where Q = 2.87 • 10 _6 s _1 is the equatorial angular velocity of the Sun at the surface, 
x = t/Rq, <f)(x,xo) = 0.5 [1 + tanh [100(x — xo)]], xq = 0.71. The distribution of the 
Coriolis number, turbulent diffusivity and the RMS convection velocity are shown in FigjTJo. 
The radial profile of the a-effect is shown in FigjlJ. 

At the bottom of the integration domain we apply the perfect conductor ( "closed") 
boundary conditions: S$ = 0, A = 0. The boundary conditions at the top are defined as the 
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Fig. 1. — Internal parameters of the solar convection zone: a) the contours of the constant 
angular velocity are plotted for the levels (0.75 — 1.05)Qo with a step of 0.025f2 , see, Eq.([7|; 
b) turnover convection time r c , turbulent diffusivity t]t, RMS convective velocity U' c ; c) the 

. The distance is measured in 



radial profile of the dynamo a-effect, a p = f]TGfiQ , see Eq. 
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following. Bearing in mind the idea of the partial escape of the toroidal flux from the Sun 
discussed in Introduction, we explore a combination of the "open" and "closed" boundary 
conditions at the top, controlled by a parameter 5. For the toroidal field we use condition: 



S^B 



[l-6)£g 



0. 



This is similar to the boundary condition discussed by Kitchatinov et al. (2000). For the 



poloidal field we apply a combination of the local condition A = and condition of smooth 
transition from the internal poloidal field to the external potential (vacuum) field: 



dA 
dr 



QAivac) 



dr 



+ (1-S)A = 0, 



(9) 



where the external potential field is : 

A^ (r, ^ = J2 a n Vl^V 2 Pn (aO 



(10) 



P^ (/x) is the associated Legendre polynomial of degree n. For the numerical implementation 
of Eq.Q, we take a one-side finite difference approximation for the radial derivative at an 
angular mesh point pif 



dA, 



dr 



3Anj — AA^-ij + An-2; 
2h r 



where h r is the radial discretization interval, and consider the expansion (10) at 
the top boundary r = r e : A^j = a n yl — /u|-P^ (a*j)- Then, we define matrices 

— iJ^Pn (fij) and Mj n = nJl — fJ%P* (/ij), with jij being collocation points of 



M. 



(a) 



nj 



P^ (see, Boyd , 2001). This procedure allows us to express coefficients a n in (10) via the 



values of potential A at the grid-points: a n 
solving it we get : 



(o)-l 



71,7 



Anj. Substituting this in Eq.(9) and 



A 



Nj = 5(3E jk + 2h r .M jk M^- 1 



{AA 



N-l jk 



-A 



N-2k) j 



where E is a unit diagonal matrix. 
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3. Results and discussion 

Parameter 5 in the boundary conditions describes a transition between the "closed'' 
(5 = 0) to "open" (5 = 1) boundaries. Physically, it controls penetration of the 
dynamo-generated fields into the outer atmosphere. 



Decreasing 8 in Eqs.(8|9) results in stronger tangential and weaker radial large-scale 



magnetic fields at the surface. While the strong toroidal magnetic field is a desired feature 
of the model, the weak radial magnetic field decreases the efficiency of the radial subsurface 
shear to produce large-scale toroidal magnetic fields. In fact, the simulations reveal that the 
critical dynamo number, C a , is greater when the penetration parameter 5 is smaller. For 
this reason we consider the case of a small deviation from the vacuum ( "open") boundary 
conditions. Moreover, in order to match the dynamo period to the solar cycle we choose 
the magnetic diffusivity parameter C v = 0.05, which is significantly lower then the value 
predicted by the mixing length theory. To demonstrate the effect of the new boundary 
conditions with the field penetration we show for comparison in Figure 2 and 3 the results 
of two runs for 5 = 0.95 (corresponding to a partial penetration of toroidal field) and 5 = 1 
(the vacuum boundary conditions). 

These results show that allowing the large-scale toroidal magnetic field to penetrate in 
to the surface layers of the Sun changes the direction of the latitudinal migration of the 
toroidal field activity and produces the magnetic butterfly diagram in a good qualitative 
agreement with solar-cycle observations. The dynamo-wave penetrates to the surface and 
propagates along the iso-surface of angular velocity in the subsurface shear layer. This is in 



agreement with the Yoshimura rule (Yoshimura, 1975) 



Dikpati et al. (2002) explored generation of toroidal magnetic fields by the fi-effect in 
the sub-surface shearlayer in the Babcok-Leighton-type dynamo models. They found that 
the phase relation between the sub-surface toroidal magnetic field and the surface radial 
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Fig. 2. — The case of 5 = 0.95 (the top boundary conditions with a partial penetration of 
toroidal magnetic fields into the outer layers of the Sun). The top panel shows the near- 
surface toroidal component of large-scale magnetic fields (contour lines) and the surface 
radial component of the field (color background). The bottom panel shows snapshots of the 
poloidal (contour lines) and toroidal magnetic field components for a half of the magnetic 
cycle. The maximum strength of the toroidal field is about 1KG. Time is in years. 




Fig. 3.- 



The same as in on the Figj2]for the vacuum boundary conditions, 5 = 1. 
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magnetic field is inconsistent with observations. We believe that this inconsistency was 
because in their model the source of the surface poloidal magnetic fields was related with 
the bottom of convection zone. Therefore, the sub-surface toroidal field that is generated 
in the sub-surface shear layer does not contribute directly to the generation of the poloidal 
magnetic field. 

Both of our simulation runs were started with the initial magnetic fields with of 
the equally mixed symmetrical and antisymmetrical realtive to the equator components. 
The evolution retains only the dipole-like parity configurations in both cases though the 
relaxation time in the penetration case of 5 = 0.95 is much longer than in the case of 
the pure vacuum boundary conditions. If we relax the confinement of the alpha-effect in 
latitude, i.e., f(9) = 1 instead of Eq.([6]), the general patterns of Figj2]are hold except that 
the maximum of the toroidal magnetic field is shifted to higher latitude w 40°. Therefore 
we can conclude that the aQ- dynamo model with the boundary conditions that allow a 
small partial penetration of the toroidal field into the outer layers of the Sun, can robustly 
reproduce the solar-cycle butterfly diagram for the near-surface large-scale magnetic field 
evolution. These results demonstrate the importance of the subsurface rotational shear 
layer in the solar dynamo mechanism. 
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4.1. Appendix 

Here, we give the definitions of the functions which were used in the model. The 
given functions describe the efficiency of the Coriolis force and the mean magnetic field 
to act on the stratified turbulence and to produce the dynamo a-effect, anisotropy of 
magnetic diffusion, turbulent magnetic pumping, magnetic quenching of the turbulent 
effects, etc. These effects are discussed in (Pipin, 2008). Functions /i^io depend on the 
Coriolis number Q* = 2t c Q q ; functions describe magnetic quenching and depend on 
= B/y/fiofnl 2 : 

Aa) 3 (, 2 , arctanQ* \ 

fw = -^(3(^ 2 + l)(^ 2 + --l)^^-((2e + l)^ 2 ^ 
^ = ^(((^ + 3) £ + 4, 2 -l)^ + l-3 £ ) 

= ^(16^-3-3(^-1)^). 

The parameter e measures the ratio between the turbulent energies of the kinetic and 
magnetic fluctuations e = ft/fiopu 2 , in the background turbulence (in the absence of 
the mean-fields). Note, in notation of Pipin(2008), the turbulent diffusivity and a-effect 
quenchning functions are defined as follows, ip v = 4>3 + 4>2 — 20i, and ip a = —3/4<j)^\ 
respectively. Expressions for 01,2,3 and are given in (Pipin, 2008). 
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